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Abstract 

The standing kink magnetohydrodynamic (MHD) modes in a zero-beta cylindrical com- 
pressible magnetic flux tube modelled as a twisted core surrounded by a magnetically twisted 
annulus, both embedded in a straight ambient external field is considered. The dispersion 
relation is derived and solved numerically to obtain the frequencies of the kink MHD waves. 
The main result is that the twisted magnetic annulus does affect the period ratio P1/P2 of 
the kink modes. 
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1 Introduction 



Transverse oscillations of coronal loops were first identified by Aschwanden et al. (1999) and 
Nakariakov et al. (1999) using the observations of TRACE (the Transition Region And Coronal 
Explorer). 

Verwichte et al. (2004), using the observations of TRACE, detected the multimode oscilla- 
tions for the first time. They found that two loops are oscillating in both the fundamental and 
the first-overtone standing kink modes. According to the theory of MHD waves, for uniform 
loops the ratio of the period of the fundamental to the period of the first overtone is exactly 
2, but the ratios found by Verwichte et al. (2004) are 1.81±0.25 and 1.64±0.23. However, 
these values were corrected with the improvement of the observational error bars to 1.82±0.08 
and 1.58±0.06, respectively, by Van Doorsselaere, Nakariakov & Verwichte (2007). Also Verth, 
Erdelyi & Jess (2008) added some further corrections by considering the effects of loop expansion 
and estimated a period ratio of 1.54. All these values clearly differ from 2. This may be caused 
by different factors such as the effects of curvature (see e.g. Van Doorsselaere et al. 2004), 
leakage (see De Pontieu, Martens &: Hudson 2001), density stratification in the loops (see e.g. 
Andries et al. 2005; Erdelyi & Verth 2007; Karami & Asvar 2007; Safari, Nasiri & Sobouti 2007; 
Karami, Nasiri & Amiri 2009), magnetic field expansion (see Verth & Erdelyi 2008; Ruderman, 
Verth & Erdelyi 2008; Verth, Erdelyi & Jess 2008) and magnetic twist (see e.g. Erdelyi &: Fedun 
2006; Erdelyi & Carter 2006; Karami & Barin 2009). 

Mikhalyaev & Solov'ev (2005) investigated the magnetohydrodynamic (MHD) waves in a 
double magnetic flux tube embedded in a uniform external magnetic field. The tube consists of 
a dense hot cylindrical cord surrounded by a co-axial shell. They found two slow and two fast 
magnetosonic modes can exist in the thin double tube. 

Erdelyi & Fedun (2006) studied the wave propagation in a twisted cylindrical magnetic flux 
tube embedded in an incompressible but also magnetically twisted plasma. They found that 
increasing the external magnetic twist from to 0.3 caused an increase in the normalized periods 
of sausage MHD waves approximately by 1—2%. 

Erdelyi & Carter (2006) used the model of Mikhalyaev & Solov'ev (2005) but for a fully 
magnetically twisted configuration consisting of a core, annulus and external region. They 
investigated their analysis by considering magnetic twist just in the annulus, the internal and 
external regions having straight magnetic field. Two modes of oscillations occurred in this 
configurations; surface and hybrid modes. They found that when the magnetic twist is increase 
the hybrid modes cover a wide range of phase speeds, centered around the annulus, longitudinal 
Alfven speed for the sausage modes. 

Carter & Erdelyi (2007) investigated the oscillations of a magnetic flux tube configuration 
consisting of a core, annulus and external region each with straight distinct magnetic field in 
an incompressible medium. They found that there are two surface modes arising for both the 
sausage and kink modes for the annulus-core model where the monolithic tube has solely one 
surface mode for the incompressible case. Also they showed that the existence and width of an 
annulus layer has an effect on the phase speeds and periods. 

Carter & Erdelyi (2008) used the model introduced by Erdelyi & Carter (2006) to include 
the kink modes. They found for the set of kink body modes, the twist increase the phase 
speeds of the modes. Also they showed that there are two surface modes for the twisted shell 
configuration, one due to each surface, where one mode is trapped by the inner tube, the other 
by the annulus itself. 

Ruderman (2007) studied the nonaxisymmetric oscillations of a compressible zero-beta thin 
twisted magnetic tube surrounded with the straight and homogeneous magnetic field taking the 



2 



density stratification into account. Using the asymptotic analysis he showed that the eigen- 
modes and eigenfrequencies of the kink and fluting oscillations are described by a classical 
Sturm— Liouville problem. The main result of Ruderman (2007), which also has been already 
obtained by Goossens, Hollweg & Sakurai (1992), was that the twist does not affect the kink 
mode. 

Karami & Barin (2009) studied both the oscillations and damping of standing MHD surface 
and hybrid waves in coronal loops in presence of twisted magnetic field. They considered a 
straight cylindrical incompressible flux tube with magnetic twist just in the annulus and straight 
magnetic field in the internal and external regions. They showed that both the frequencies and 
damping rates of both the kink and fluting modes increase when the twist parameter increases. 
They obtained that the period ratio P\jPi of the fundamental and first-overtone for both the 
kink and fluting surface modes are lower than 2 (for untwisted loop) in presence of the twisted 
magnetic field. 

In the present work, our aim is to investigate the effect of the twisted magnetic annulus on 
the period ratio P1/P2 of kink MHD waves in the coronal loops observed by Verwichte et al. 
(2004) deduced from the TRACE data. This paper is organized as follows. In Section 2 we use 
the asymptotic analysis obtained by Ruderman (2007) to derive the equations of motion. In 
Section 3, using the relevant boundary conditions, we obtain the dispersion relation. In Section 
4, we give numerical results. Section 5 is devoted to conclusions. 

2 Equations of motion 

The linearized MHD equations for a zero-beta plasma are 

' [(Vx5B)xB + (VxB)x5B], (1) 

V x (Sv x B), (2) 



dt 4vrp 

d5B 



dt 

where 5v and (5B are the Eulerian perturbations in the velocity and magnetic fields; p, is the 
mass density. 

The simplifying assumptions are as follows. 

• The background magnetic field is assumed to be 
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Bi = (0,A ir ,B zi (r)), r<a, 

B = (0,A r,B z0 (r)), a < r < R, (3) 
B c = (0,0, B ze ), r>R, 



where A- x , A®, B zc are constant and a, R are radii of the core and tube, respectively. 

From both the equilibrium equation, i.e. = and the continuity condition of the 

magnetic pressure across the boundaries of the tube, i.e. B? (a) = B^(a), Bq(R) = B 2 (R), 
the z-component of the equilibrium magnetic field can be obtained as 

B^r) = Bg + 4V-2r 2 ), 

BUr) = Bi + Al{a" - 2r 2 ), (4) 
Bl = B> + AUa 2 -R 2 ), 

where Bq is an integration constant. The above magnetic field configuration in the absence 
of the annulus is the same as the background magnetic field considered by Ruderman 
(2007). 



3 



• p is constant over the loop but different in the interior, annulus and exterior regions and 
denoted by p\, po and p e , respectively. 



• Tube geometry is a circular with cylindrical coordinates, (r, <f>, z). 

• There is no initial steady flow over the tube. 

• t-, (p- and z- dependence for any of the components 5v and <5B is exp {i(m</> + k z z — u>t)}. 
Where k z = lir/L, L is length of the tube, and I = (1,2, ■ ■ •), m = (0, 1,2, • • •) are the 
longitudinal and azimuthal mode numbers, respectively. 

Like Ruderman (2007), we define e = ^ <C 1 which is in good agreement with the ob- 
servations. Following the second order perturbation method in terms of e given by Ruderman 
(2007), solution of Eqs. (l)-(2) in terms of P = ^-P, the Eulerian perturbation in the magnetic 
pressure, and £ r = —5v r /iuj, the Lagrangian perturbation in the radial displacement, for the 
twisted regions yields 



where F = k z + m-^. Equations (5) and (6) are same as Eqs. (19) and (21), respectively, in 
Ruderman (2007). 

In the interior and annulus regions, solutions of Eq. (6) are 

{Q<T m ^ T <C Ob 

p r m-i + 7r -™-i ; a <r <R. ^ 

For the exterior region, we obtain 



47 r dP 

" "PBS <9) 

Equations (8) and (9) are same as Eqs. (26) and (25a), respectively, in Ruderman (2007). In 
the exterior region, r > R, the waves should be evanescent. Solutions are 

P(r) = sK m (k'r), k' 2 > 0, (10) 



Cr(r) = -e^K'Jk'r), (11) 

where K m is the modified Bessel function of the second kind and a prime on K m indicates a 
derivative with respect to its appropriate argument. The coefficients a, (3, 7 and e in Eqs. (7) 
and (10) are determined by the boundary conditions. 
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3 BOUNDARY CONDITIONS AND DISPERSION RELATION 



Following Ruderman (2007), at the perturbed tube boundary the plasma displacement in the 
radial direction and the magnetic pressure have to be continuous as 



e 



r0 
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r0 



=R 



--R' 



(12) 



B 2 



P ~ ISS^rO 



£r0 



(13) 



r=_R 



Using this boundary conditions and the solutions obtained in the previous section, the dispersion 
relation is derived as 



/a\2m„ i 4ir R K' m (k' R) „ { 
\RJ ~ ~Bjk> K m (k'R)^ r 

■4tt RK' m (k'R)„ „ 



i? 

+ 



■=•0 



a \ 2m 



(I) 



B$k'K m (k'Ry 



R 

/ a \2m 



= 0, 



(14) 



with 



B 2 k 2 
Air 



+ 



47rm 



(2% + m4 J )(l-m), 



(15) 



H m = --(p u, 2 - ^) + -^(2B k z + mi )(l + m), (16) 

m V 47T / 47T771 

where in 3^, the superscript j stands for i and corresponding to the interior and annulus 
regions, respectively. 

Note that if we remove the annulus region, i.e. setting a = R, then the four boundary 
conditions, Eqs. (12)-(13), reduce to two boundary conditions and finally the dispersion relation, 
using the thin flux tube approximation for K m (x) oc x~ m at small x, yields to 

u? = Cl{kl + ^^^B Q k z + Am)}, (17) 

B 2 

where Of: = 2 -k( p -Xp ) ' Equation (1?) xs sam e as Eq. (40) in Ruderman (2007). The main result 
of Ruderman (2007) is that the twist does not affect the kink modes and Eq. (17) shows that we 
get the same frequencies as in the case that A\ = 0. This result also has been already obtained 
by Goossens, Hollweg & Sakurai (1992). Note that Eq. (14) shows that even in the presence of 
the annulus, Aq / 0, the internal twist does not affect the kink (m = 1) modes. Because the 
internal twist, Ai, only appears in Eq. (15) and when m = 1 then it has no contribution. 

In the next section, using the numerical solution of the dispersion relation, Eq. (14), we 
show that the twisted annulus region, which adds a new boundary to the system, does affect 
the frequencies of the kink modes. 
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4 Numerical Results 



As typical parameters for a coronal loop, we assume L = 10 5 km, a/L = 0.01, p c / pi = 0.1, 
Po/Pi = 0.5, pi = 2 x 10~ 14 gr cm" 3 , B = 100 G. For such a loop one finds v A . = ^gL= = 2000 

km s _1 , uJAi '■= = 0.02 rad s _1 . 

The effect of twisted magnetic annulus on the frequencies u> is calculated by the numerical 
solution of the dispersion relation, Eq. (14). The results are displayed in Figs. 1 to 4. Figures 
1 to 2 show the frequencies of the fundamental and first-overtone I = 1, 2 kink (m = 1) surface 
modes with radial mode numbers n = 1, 2 versus the twist parameter of the annulus, B^/B z := 
4^, and for different relative core width a/R = (0.65,0.9,0.99). Figures 1 to 2 reveal that: i) 
for a given a/R, the frequencies increase when the twist parameter of annulus increases. The 
result is in good agreement with that obtained by Carter & Erdelyi (2008) and Karami & Barin 
(2009). Note that the existence of the two surface modes labelled by (n = 1,2), corresponding 
to the two boundaries located at r = a, R, are also in accordance with Carter & Erdelyi (2008) 
and Karami & Barin (2009). ii) For a given n and a/R, when the longitudinal mode number, 
I, increases, the frequencies increase, iii) For a given a/R and B$/B z , when the radial mode 
number, n, increases, the frequencies increase, iv) For n = 1, when a/R goes to unity then the 
frequencies become independent of B$/B z . Therefore in the absence of the annulus, the twist 
does not affect the kink modes. This is in good agreement with those obtained by Goossens, 
Hollweg & Sakurai (1992) and Ruderman (2007). v) For n = 2, when a/R goes to unity exactly 
then the frequencies are removed. This is expected to be occurred because for a/R = 1 we have 
only one surface boundary corresponding to one surface mode. 

The period ratio P1/P2 of the fundamental and first-overtone, I = 1,2 modes of the kink 
(m = 1) surface waves with n = 1, 2 versus the twist parameter of the annulus is plotted in Figs. 
3 to 4. Figures 3 to 4 show that: i) the period ratio P1/P2 with increasing the twist parameter 
of the annulus, for n = 1 decreases from 2 (for untwisted loop) down to a minimum and then 
increases. Whereas for n = 2, it decreases from 2 and approaches below 1.6 for a/R = 0.5, for 
instance. Note that when the twist is zero, the diagrams of P1/P2 do not start exactly from 
unity. This may be caused by the radial structuring (po / pi, p e / pi). But for the selected 
thin tube with a/L = 0.01, this departure is very small, O(10 -4 ), and doesn't show itself in the 
diagrams (see McEwan et al. 2006). ii) For a given B^/B;,, the period ratio P1/P2 for n = 1 
increases and for n = 2 decreases when the relative core width increases. Figure 3 clears that for 
kink surface modes (m = l,n = 1) with a/R = 0.5, for both B^/ B z =0.0ll and 0.015 the ratio 
P1/P2 is 1.82. This is in good agreement with the period ratio observed by Van Doorsslaere, 
Nakariakov & Verwichte (2007), 1.82±0.08, deduced from the observations of TRACE. See also 
McEwan, Diaz & Roberts (2008). 

5 Conclusions 

Oscillations of standing kink MHD surface waves in coronal loops in the presence of the twisted 
magnetic annulus is studied. To do this, a typical coronal loop is considered as a straight 
cylindrical compressible zero-beta flux tube with magnetic twist in the internal and the annulus 
and straight magnetic field in the external region. Using the perturbation method given by 
Ruderman (2007), the dispersion relation is obtained and solved numerically for obtaining the 
frequencies of the kink modes. Our numerical results show that 

i) for a given relative core width, frequencies of the fundamental and first-overtone / = 
1,2 kink (m = 1) surface modes with radial mode numbers n = 1,2 increase when the twist 
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parameter of the annulus increases; 

ii) when the relative core width, a/R, goes to unity exactly then the kink (m = 1) surface 
modes with n = 1 become independent of the twisted annulus and the second modes labelled 
by n = 2 are removed from the system; 

iii) the period ratio P1/P2, for the kink [m = 1) surface modes with n = 1,2 is lower than 
2 (for untwisted loop) in the presence of the twisted magnetic annulus. The result of P1/P2 for 
the kink (m = 1) surface modes with n = 1 is in accordance with the TRACE observations. 
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Kink Surface Modes (n=l,m=l,l=l) 
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Figure 1: Frequencies of the fundamental and its first-overtone kink (m = 1) surface modes 
with radial mode number n = 1 versus the twist parameter of the annulus, B^,jB z = ^§^, for 
different relative core width a/R = 0.65 (dash-dotted), 0.9 (dashed) and 0.99 (solid). The loop 
parameters are: L = 10 5 km, a/L = 0.01, p c / p\ = 0.1, po/p\ = 0.5, p\ = 2 x 10 -14 gr cm~ 3 , 
Bq = 100 G. Frequencies are in units of the interior Alfven frequency, wa; = 0.02 rad s -1 . 



Kink Surface Modes (n=2,m=l,l=l) 




Kink Surface Modes (n=2,m=l,l=2) 
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Figure 2: Same as Fig. 1, for the kink (m = 1) surface modes with radial mode number n = 2. 
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Kink Surfase Modes (m=l,n=l) 
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Figure 3: The period ratio Pi /Pi of the fundamental and its first-overtone kink (m = 1) 
surface modes with radial mode number n = 1 versus the twist parameter of the annulus for 
different relative core width a/R = 0.5 (dotted), 0.65 (dash-dotted), 0.9 (dashed) and 0.99 
(solid). Auxiliary parameters as in Fig. 1. 
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Figure 4: Same as Fig. 3, for the kink (m = 1) surface modes with radial mode number n = 2. 
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